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Abstract. We present a two-layer accretion disc model developed to simultaneously fit optical long baseline vis- 
ibilities and spectral energy distributions of T Tauri accretion discs. This model allows us to access easily the 
physical conditions in the disc as the mid-plane or the surface temperature. 

Our model includes viscous heating, absorption of stellar irradiation, and thermalisation with the surrounding 
medium. The disc is modelled with concentric cylinders for which the vertical radiation transfer is computed using 
two layers with vertically averaged temperatures: the outer layer is heated by the stellar irradiation and by the 
inner layer, and the inner layer by viscous dissipation and by the outer layer. We investigate three prescriptions for 
the geometrical thickness of the disc: it is either proportional the scale height (model 1), given ad hoc (model 2), 
or zero (model 3). We then derive the disc structure in the case of the a and (3 viscosity prescriptions, as well as 
for various optical thickness regimes of the disc. 

This analytical model allows us to disentangle regions where the mid-plane temperature and the effective tem- 
perature are dominated by accretion from regions dominated by reprocessing of stellar light. In the case of 
a-prescription, we find that the structure of model 2 gives predictions very close to those of numerical simulations 
from previous authors. 

From the disc structure, we derive the spectral energy distributions, images and interferometric visibilities. We 
analyse the influence of the disc parameters on the resulting structure and on the observable outputs. We apply 
our model to interpret consistently the spectral energy distributions and visibilities of SU Aur and FU Ori for 
which interferometric data are available, and that are not known to be part of a multiple system. We were not 
able to derive a consistent fit for T Tau North, which might come from caveats in the flux correction from its 
South component, but were able to separately derive fits for its spectrum and its visibilities. 

We find that even a single interferometric measurement at one infrared wavelength can bring a very strong 
constraint on disc models. We predict that future massive interferometric observations of accretion discs will 
provide a breakthrough in the understanding of accretion disc physics. 

Key words. Stars: pre-main-sequence; circumstellar matter - Accretion; accretion disks - Radiative transfer - 
Methods: analytical; numerical; data analysis 



1. Introduction 

Since the initial m o dels o f visc ous accretion discs by 
IShakura fc SunvaevI l|l973|) and ILvnden-Bell fc Pringld 
the physics of the close environment of T Tauri 
stars (TTS) has been extensively studied in order to in- 
terpret their spectral energy distribution (SED). For the 
sake of simplicity, models traditionally separated discs into 
two categories, sometimes called active discs, on one hand, 
in which viscous dissipation is predominant, and passive 
discs, on the other hand, for which irradiation by the 
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central star is the main heating process. Early models 
used quasi-Keplerian steady accretion discs, assumed to 
be geometrically thin for a wide range of accretion rate; 
they predicted a fixed slope for the infrared spectrum: 
A F\ ~ A" 4/3 . However many TTS present flatter SEDs, 
and disc flaring was among the first attempts to explore 
disc ve rtical structure as an explanation for such S EP flat- 
tening l|Adams et al.ll987UKenvon fc Hartmann1ll987tk in 
a flared disc, the surface of the remote parts is tilted to- 
ward the star and gets more stellar light than forecast by 
the standard model, resulting in a warmer disc further 
away from the star. Since then, several models have been 
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proposed in order to explain both standard SEDs and flat- 
ter ones. 

An analytical study of the radiative transfe r in the 
vertic al structure of discs was fi rst carried out bvlHubenyl 
l|l990(l for active discs, then bv lMalbet fc BertoutI l|l99lL 
hereafter paper I) for passive discs. In the latter model, 
the topmost layers of the disc, illuminated by the star, are 
hotter than the disc photosphere, re sulting in excess con- 
tinuum and line emission. Later on, IChiang fc Ooldreicbl 
l)l997lll999]) used a simplified two-layer passive disc model 
based on the same super-heating mechanism and derived 
SEDs, confirming conclusions of paper I and produc- 
ing results consisten t with observations. More recently, 
iMalbet et alJ l|200li hereafter paper II) generalised the 
analytical study of paper I to discs heated by several pro- 
cesses, and used its formalism to derive the the vertical 
structure of active discs. 

On the other hand, numerical integration of the equa- 
tions of radiative transfer was carried out by various au- 
thors in order to derive the vertical structur e of accre- 
tion d is cs with fewer a p riori approximations. Bell fc Lin 
lll994lk iBelT et al. | lll997l) developed an a ctive disc model 
in order to e xplain FU Orionis outbursts: iD'Alessio et alJ 
l)l998l Il999|) dealt with the more general case of a disc 
heated both by viscous heating and stellar irradiation. 

From a general point of view, all these studies pre- 
dict spectra consistent with observations, but they have 
rarely been checked consistently against the spatial in- 
formation revealed by recent optical and infrared high 
angular resolution imaging. Recently, the advent of op- 
tical in terferometry has se t newer constr a ints on dis c 
models. IMalbet et al.l fjflflft: lAkeson et all lElOfl Eool: 
Mal bet fc Be rgcr ( 2002) obtained the first visibility mea- 
surements of TTS and FU Orionis stars. However, they 
failed to consistently fit both SEDs and visibilities with a 
standard disc model: most of the time, the disc parameters 
derived from the SED data are in disagreement with those 
derived from visibility data. No other attempt to compare 
self-consistent disc models and interferometric measure- 
ments has been carried out so far for low mass pre-main 
sequence stars. 

In this paper, we tackle the issue of analytically de- 
scribing a disc in presence of the two main heating pro- 
cesses, viscous heating and stellar irradiation, the latter 
requiring a correct description of the flaring. This model 
suits the TTS and FU Ori-type stars for which viscous 
heating cannot be ignored; for more massive stars (H erbig 
Ae/Be) it has been shown i Dullemond et al1l200lh that 
viscosity, as a heating mechanism, can be ignored, so that 
our model is not relevant. In Sect. [2 we present a two- 
layer version of the model developed in paper I and carry 
out an analytical determination of the structure of the 
disc. We derive a set of equations giving the mid-plane 
temperature and the flaring index, from which the whole 
structure can be determined. In Sect. [31 we briefly present 
the numerical approach. In Sect.jTj we compare the results 
of the model with other models and analyse the influence 
of some disc parameters on the observables. In Sect.0 we 



apply our model to the few low-mass young stellar objects 
(YSOs) observed in optical interferometry. 

2. Model description 

2.1. The two-layer structure 

The st andard accretion disc model bv lShakura fc Sunvaevl 
determines the disc emission with its effective tem- 
perature and does not take the vertical temperature profile 
into account. However, this approximation fails when the 
disc surface is super-heated by stellar irradiation (paper 
I): the illumination by the star is predominant in a opti- 
cally thin outer layer of the disc while the optically thick 
part of the disc is governed by the radiative transfer of 
thermal radiation. This phenomenon results in emission 
at different temperatures and has a strong incidence on 
the SED. 

To further explore the resulting structure of the disc, 
we use a two-la yer vertical structure. This approach was 
first proposed bv lChiang fc Goldreichl l)l997|) in the case of 
passive discs, and we apply it here in a more general case 
where viscous dissipation, stellar irradiation, and thermal- 
ization with the surrounding medium are taken into ac- 
count. In the case of a passive disc the disc physical prop- 
erties varia tions with distance to the s tar follow a unique 
power law (|Chiang fc OoldreichllT997l) . In an active disc, 
where the temperature is dominated by viscous dissipa- 
tion, this is no longer the case (paper II, see Fig. 3), due 
to the sensitivity of the vertical structure to the disc ma- 
terial opacity. 

In this paper, we present a model of disc that can be 
seen as a simplification of the formalism of paper II. In 
this model, we consider an optically thick disc heated both 
by stellar irradiation and viscous dissipation. An optically 
and geometrically thin outer layer is directly heated by the 
star and by the inner layer; the optically thick inner layer 
is heated by viscous dissipation and by the outer layer. In 
the outer layer we use a vertically averaged temperature to 
solve the radiative transfer, and in the inner layer we use 
the mid-plane temperature. A radial slice of such a disc is 
represented in Fig. ^ The inner and outer layers have the 
optical thicknesses r; and r , and the temperatures 7] and 
T e respectively. 

The incidence of stellar radiation onto the disc is given 
by the angle cj>, which is an average of this incidence over 
the surface of the star. The height of the outer layer is Ti. 

In this section, we focus our attention on two quanti- 
ties: the mid-plane temperature and the flaring angle. The 
structure of the disc can be derived from these two quanti- 
ties; for instance, the mid-plane temperature governs the 
scale height and the flaring angle the surface tempera- 
ture. In order to obtain 7] and 7 we proceed in four steps: 
in Sect. 12.21 we derive the effective temperatures T v , T T 
and T a corresponding to viscous heating, stellar light re- 
processing and thermalization with the disc surrounding 
medium. In Sect. EH equations of transfer are obtained, 
giving the temperatures 7] and T c as a function of the 
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Fig. 1. Radial slice of a two-layer disc. 

optical depths and t c and the effective temperatures 
T V) T r and T a . In Sect. 12.41 the structure of the disc (scale 
height, column density, optical depth, etc.) is connected to 
the temperatures XI and T c . In Sect. 12.51 the results from 
the first three steps are combined to derive the mid-plane 
temperature 7] and the flaring (variation of TL/r with r) 
in a set of coupled equations. 



as in a stellar atmosphere within the Eddington approxi- 
mation. Concerning the temperature in the mid-plane of 
the disc, as showed in paper II, the stellar atmosphere for- 
malism cannot be used as is. A corrective term At; must 
be applied to the optical thickness ti: 



An). 



(5b) 



For the sake of simplicity we ignored the usual constant 
term 2/3 before the large optical depth. 

If reprocessing only occurs, the outer layer is super- 
heated because it is optically thin: the stellar irradia- 
tion is dissipated along a slanted path of optical depth 
unity in the visible, i.e. much less than one for the re- 
pro cessed IR radiation in the vertical direction. We fol- 
low IChiane&GcJdrcic hl lll997l) and state that both layers 
present a vertically averaged temperature. The outer layer 
catches the flux T 4 and emits half of it upward, while be- 
ing optically thin, so that 



2.2. Heating sources 
iLvnden-Bell fc Pringlel 



showed that the effective 
temperature of a geometrically thin active disc is 



(77T 



where / is given by 



/ = 1 



(1) 



(2) 



M , J7k and are respectively the accretion rate, the 
Keplerian angular velocity at radius r, and the radius of 
the star. 

We suppose that the surface of the disc presents a 
mean albedo A constant over the whole disc. The effective 
temperature associated with stellar heating then reads 



rp4 



0T C 



cq ' 



(3) 



T eq is the Local Thermal Equilibrium temperature of a 
body with an albedo A at a distance r of the star. It is 
given by 

T c \-(l-A)^ (^) , (4) 
where T* is the effective temperature of the star. 



2.3. Equations of transfer 

Since the layers are assumed to be isothermal, the radia- 
tive transfer has a simple form, depending only on the 
optical depths of the layers and their temperatures. 

If viscous dissipation only occurs, we expect a surface 
temperature of 



rpl 



T 



(5a) 



rp4 



2r c 



(6a) 



and emits the other half down to the optically thick inner 
layer, hence 



T: 4 = ±1-. 



(6b) 



In the case of thermalization with the disc surrounding 
medium, the temperatures are 



T? = T 4 , 

rpi rp4 



(7a) 
(7b) 



In the general case, the three heating contributions 
considered above occur together. The superposition prin- 
ciple can be applied to the specific radiative intensities; 
that is the radiation at a given location of the disc is the 
sum of the radiations corresponding to different heating 
processes (see paper II). Summing over Eqs. (0 El 13 we 
derive 



rp4 I , 



rp4 rp4 



2 

_ r 

2r n 



-T 4 (r 

- 1 v 

2 ' 



ArO, 



(8a) 
(8b) 



This result can also be obtained from the equations of 
transfer 28-31 of paper II using the following approxima- 
tions: (i) the surface is optically thin, (ii) both layers are 
isothermal, (iii) viscosity only occurs in the inner layer, 
(iv) diffusion is ignored, (v) no stellar radiation penetrates 
the inner layer, and (vi) the vertical optical thickness of 
the disc is larger than unity. 
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2.4. Structure of the disc 

If we ignore self-gravity, the vertical gravitational field is 
proportional to the distance z from the mid-plane. Since 
the inner layer is isothermal, the density reads 



(9) 



p(z) = piexp ( -^-j 



where the density scale height hi is linked to the sound 
speed c; and the gravitational field g z by 



hi = C; 



dg z 
dz 



-1/2 



(10) 



We consider the gravitational temperature T g at which 
the thermal energy of a particle balances the gravitational 
energy, given by 



GM»m H 



where mjj is the atomic mass of hydrogen. Then, 



(11) 



(12) 



where pi is the mean molecular mass in the mid-plane of 
the disc. 

In the standard disc model by iLvnden-Bell fc Pringlel 
the mass column £ is linked to the uniform accre- 
tion rate M by 



M 
Sirvi 



f, 



(13) 



where Vi is the vertically averaged viscosity. 

In order to go f urther, we need a prescrip tion for the 
kinematic viscosity. IShakura & Sunvaevl ( 1973h use the so- 
called a-prescription: 



(14) 



where a is a dimensionless factor ranging from 10 -4 to 
10 _1 . If we ignore the thin outer layer, the viscosity is 
uniform over the entire height because of the isothermal 
assumption. So, Eq. (|13[1 . together with Eq. Q14JI. becomes 



Ti 



The mass column also meets 

£ = 2 / p(z)dz. 
Jo 

so that, with Eqs. flUI^I 



Pi = Po 



Ti 



, Po 



2 Mf 



7T 6^x0! 7rr 3 



(15) 



(16) 



(17) 



In order to keep an analytical description, we assume 
that the opacity k locally follows a power law: 



Irj-rni 



K(p, T) = kp l T 



(18) 




Fig. 2. Incidence angle of stellar radiation onto the disc 

We now assume that the outer layer has an optical 
depth much smaller than that of the inner layer, that is 
t <C Ti . This is indeed true in optically thick regions of the 
disc, where Ti 3> 1 and T e <C 1. In optically thin regions, 
the validity of this approximation must be self-consistently 
checked after computation. With such an approximation, 
the optical depth of the disc can be written as if the outer 
layer were absent: 



^.+00 

Ti = / p(z)k(z) dz, 
Jo 

so that 



Ti = 



-kp^T^h 



2(1 + 1 



(19) 



(20) 



The corrective term At; depends on the distribution $ of 
the viscous heating. Paper I gives 



j7°°p(CMC)dC 



r+00 

At; = / p(z)K{z)d{z)dz. 
Jo 

After some calculations we finally obtain 



^ - At; = 



with 



1 



-WT^hi 



1 



VTTT vT+2 



(21) 
(22) 

(23) 
(24) 



We call u> the ratio of opacity for radiation reprocessed 
by the outer layer to opacity for stellar radiation, uj is 
smaller than unity. With such a definition, the optical 
depth of the outer layer is 



LU(j). 



(25) 



The incidence angle of stellar radiation onto the disc 
is related to the flaring of the disc. Yet, in a flat disc, this 
incidence is non-zero because of the vertical extent of the 
star. FigureElshows that <fi is the sum of two contributions: 
the mean extent of the star 4>* and the tilt of the disc 
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surface toward the star <fif. Since the star is not a point- 
like source, we need to integrate over the stellar surface to 
find We introduce the flaring angle as the logarithmic 
variation of the disc thickness H/r with respect to the 
radial location: 



ith 



_ dlnH _ 

1 ~ ahT7 ~ ' 



(26) 



7 represents the variation of the opening angle with the 
distance to star, i.e. H/r ~ r 1 . If the opening angle of the 
disc is constant, e.g. in a flat disc, 7 is zero. With such a 
definition, the incidence angle is 



k + . 



with 



H 

1—, 

r 

4 r* 
37T r ' 



(27a) 

(27b) 
(27c) 



The expre ssion of 4>* is an asymptot ic value at large r, 
as used in lChiang fc Goldreic Then, we need to 

connect H to the previously determined scale height hi 
and use the closure relation 



H = Qh h 



(28) 



IChiang fc Goldreichl (Il997l) notice that Q is a slowly vary- 
ing quantity close to 4 and imposed Q = 4 over the disc. 
iD'Alessio et alJ (^999) consistently determine the location 
of the outer layer, and we used their results for the fiducial 
model to derive that Q ranges from 3 to 5 (see Fig. I5h.|l . 
Finally, 



(29) 



2.5. Implicit equations governing the structure 

First, we use the previous developments to derive the mid- 
plane temperature 7] as a function of the disc parameters. 
Second, we derive the flaring index on which the temper- 
ature also depends. Last, we determine a boundary condi- 
tion on H that allows to derive it over the while disc using 
H/r oc r 7 . 




2.5.1. Mid-plane temperature 



All quantities of Eq. 18a|l have been expressed in terms of 
Zj and /Ltj, of the physical constants of the disc, and of the 
quantities 7 and Q. Therefore, we derive an equation on 



T 4 = 5>(r,TKr)) 4 



h(r,Ti) 4 

t 3 (r,Ti) 4 = ~T? 
; 4 (r,Ti) 4 =T a 4 . 



3 1 _4_l+|i„ 1 -fZ+m.-l 



1. 



(30b) 

(30c) 

(30d) 
(30e) 



Equation H30a(l states that the mid-plane temperature 
is the sum of different heating processes, tk represents the 
thermal contribution of the fc-th process. t\ (r) is the mid- 
plane temperature for an optically thick disc dominated by 
viscous heating; t<z(r) is the mid-plane temperature for a 
flared disc dominated by stellar irradiation, if the geomet- 
rical extent of the star is ignored due to the distance; t^(r) 
is the mid-plane temperature for a flat disc dominated by 
stellar irradiation; £4 (r) is the mid-plane temperature in 
a disc thermalised with the ambient medium. The radial 
and thermal dependencies of the contributions tk are 



h(r,Tif ~ r - 31 - 9 / 2 /(r) ,+2 Tf 3/2,+m - 1 , 

t2(r,T{)*~r- a ' 2 TV 2 , 
t 3 (r,Ti)^r- 3 , 
i 4 (r,Ti) 4 ~ 1. 

2.5.2. Flaring index 



(31a) 
(31b) 
(31c) 
(31d) 



The case of the flar i ng is stlightly more complicated. In 
IChiang fc Goldreichl l|l997h . the disc thickness was both a 
power law of the distance to the star (H oc r 9 / 7 at large 
radii) and proportional to the scale height hi. In presence 
of viscous dissipation, it is no longer possible to ensure 
that both properties are met, since the scale height is not 
a simple power law. We investigate three possibilities: 

1. The thickness is proportional to the scale height 
(model 1). 

2. The thickness is given by the same law as in the passive 
disc (model 2). 

3. The disc is flat, which occurs if the dust grains settle 
in the mid-plane (model 3). 

The flaring index of model 1 can be expressed as the 
average of the specific flaring indices 7^ imposed by the 
various heating processes weighted by their characteristic 
temperatures tk (see Appendix iBl for a demonstration). 



(30a) 



(32a) 
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with 



obtain 



71 = 

72 = 

73 = 

74 = 



1 + 2 



1 - 31 - 2m 

20 + 61 - Am + 20 + 61 - Am 
2 

r 
i 



a- 



(32b) 
(32c) 



where 

1-/ 



9 = 



f 



(32e) 



(33) 



arises from the departure of viscous heating from the ra- 
dial power law in r -3 . This effect occurs close to the star 
and is characterised by / (Eqs.njEJ- 7i is the flaring index 
for an optically thick active disc, therefore depending on 
the opacity law; 72 the flaring index in a passive disc, if 
the geometrical extent of the star is ignored; 73 the flaring 
index in a flat passive disc; 74 the flaring index in a disc 
therm alised with the ambie n t me dium. 

In IChia.ng fc Qoldreichl <|1997|) . the thickness meets 
TL oc r 7 , with 7 = 1/8 close to the star and 7 = 2/7 
at large radii, depending whether the dominant effect in 
stellar light reprocessing is the extent of the star (£3) or 
the disc flaring (£2). We have assumed in model 2 that, 
even in the presence of viscosity, we can still write 



7 



724 



73t| 



^9 + 



(34) 



which means the disc thickness TL is only governed by the 
irradiation terms (£2 and £3) and varies from 1/8 to 2/7 
from the central parts to the outer ones. 
In model 3, the disc is flat, that is 7 = 0. 

2.5.3. Disc thickness on the outer edge 

In mo del 1 we have assumed, like in lChiang fc Goldreichl 
lll997ft . that the closure relation TL = Qh\ applies, with Q 
constant over the disc. The authors took Q — 4, which is 
seemingly a good approximation for low-mass discs (M « 
lO _8 M0/yr) but which might no longer be correct over 
the wide range of masses that we shall consider (M rj 
10" 8 -10- 4 M /yr). For this reason we prefer to determine 
the value of Q, and chose to carry out the calculation at 
the outer edge of the disc. 

In model 2, the thickness is no longer proportional to 
the scale height, so that we have to determine TL with the 
relation ffi^r = 7+1, that requires a boundary condition, 
and therefore determine TL at the other edge of the disc. 
We perform this calculation by computing Q = TL/h^. 

So, we need to determine the value of Q at the outer 
boundary of the disc in both models, which leads us to 
approximately solving the vertical hydrostatic equilibrium 
and the radiative transfer of the incoming radiation. We 



\ S„{Q) 



2 h (w/^ (Q) , 



jQuj V T\ 



(35) 



(32d) where 5h and S p are corrections due to the deviation of the 
gravity field from the linear law g z cx z when the elevation 
becomes of the order of r: 



-I/a' 



S P (Q) = 2 f-j [l - (1 + {Qh/rf) 
Sh(Q) = ((1 + (Qhi/r 2 )) 3/2 . 
A demonstration can be found in Appendix El 



(36) 
(37) 



3. Numerical approach 

Though analytical, the expressions of T t and 7 do not 
allow direct determination, because the equations giving 
7 and T; are coupled and the opacities do not follow a 
power law over the whole disc. This section first describes 
the numerical method we used for the derivation of the 
structure. Then, we briefly explain who the observables 
(SEDs and visibilities) have been derived. 

3.1. Opacities 

Because the heating and the flaring of the disc is depen- 
dent on opacities, their adequate description is manda- 
tory. The computation of the SED and image requires 
monochromatic opacities that we took from two sources: 

— At low temperatures (T < 1800 K) the opacity is dom- 
inated by the dust. We use opacities computed by 
Hcn ning fc Stognienkol l)l996(l for a mixture of inhomo- 
geneous aggregates; they include continuum emission 
and the silicate feature in the mid-infrared. 

— At higher temperatures (T > 1800 K) the opacity 
comes from to the gas. We could not find such opaci- 
ties; instead, we used piecewise power laws K\(p, T) = 
Ko(p,T)X~ p , and chose (p,T) and q, so th at they ap- 
proxim ately match Rosseland opacities bvlBell fc Lin 
( 1994) and Planck opacities bv lHenning fc Stognienkol 
(1996). Despite of this poor description, the influence 
on the observables is small because it only concerns 
the inner disc, where both optical depth (r 3> 1) and 
thickness (TL <§C r) are not critical. 

These opacities were smoothed in order to avoid numerical 
instabilities between the different domains; it is justified 
since we use a vertically averaged structure in both layers. 

The structure requires the Planck and Rosseland opac- 
ities, as well as the determination of the opacity ratio 
lu. They have been determined from the monochromatic 
opacities as explained in Appendix ^ In order to keep 
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the power law opacities formalism, we determined the lo- 
cal equivalent kr oc p T m with 



dlog «R 
<91ogp 

dlog k r 
dlogT ■ 



(38a) 
(38b) 



3.2. Physical conditions in the disc 

The disc is divided in concentric cylinders and the two- 
layer formalism is applied for each one: it consists in com- 
puting tfc, refining the values for 7] and 7 (Eqs. I30I32|) . 
computing the corresponding opacities with the associ- 
ated / and m (Eq. I38|) . deriving TL, and iterating until 
convergence. In model 1, the computation is carried out 
in ascending order of radius Cylinders shadowed by in- 
ner ones do not receive stellar light and are determined 
with the same equations, but using T* = OK. H. is de- 
termined as H = Qh with Q constant; the self-consistent 
determination of Q at the outer edge has not been imple- 
mented so far. In model 2, Q is determined at the outer 
edge (Eq. 128 jl . then the computation is performed from 
outer to inner cylinders. For each cylinder, TL is given by 
the thickness of the enclosing cylinder Tt out using 



H = H° 



(—) 

y j,out J 



1+7 



(39) 



The validity of the derived structure was checked with 
two tests. On one hand, when setting 7 = and t = 1, the 
results are co nsistent with those of the standard disc as 
determined in lBertout et al.l 0988). On the other hand, 
the model locally determines the flaring index (Eq. I32fl 
without using the radial structure, and we checked that 
the derived radial structure fulfilled ^P- = 7. 



3.3. Spectral energy distributions 

In order to compute a synthetic image or a SED, we di- 
vide each disc annulus in angular sectors. For each sector 
we determine the emergent flux in the observer's direc- 
tion as the sum of three contributions: the contributions 
of the outer layers F% , of the inner layer F x , and of sim- 
ple isotropic scattering of stellar radiation by the surface 
of the disc F|. These fluxes depend on the wavelength- 
dependent optical thicknesses t x and r| of the inner and 
outer layers. The incidence of the line of sight onto the 
visible surface of the disc is given by its angle i, and its 
incidence onto the opposite surface by i' . 



ith 



n = 



Ft 



Asm(cj)) (t 



(7) ^ (T * 



(40b) 

Bx(T e ), 

(40c) 

(40d) 



where B\ is the black-body function. T? ot (i) is the tem- 
perature of the inner layer at optical depth cos i at wave- 
length A if t x > cosi, and its mid-plane temperature in 
the other case. 

This determination allows a fast computation but is 
only an approximation. It is only valid when the outer 
parts of the disc shadows neither the star nor other parts 
of the disc; therefore a more adequate radiative transfer 
determination must be carried out for edge-on discs, which 
is beyond the scope of this paper. 

As a test, we have checked that the bolometric flux 
predicted by the SED is consistent with the one predicted 
by the effective temperatures of the heating processes. We 
noticed discrepancies of 2-3% in most cases, up to 8% in 
a few ones; they are connected to the use of the mean 
opacities kr and up in the structure calculation. 

3.4. Visibilities 

The image and visibilities are determined using the same 
method as for the SED. However, optical visibilities at 
long baselines can be both sensitive to large-scale struc- 
tures up to 1" (e.g. Palomar Testbed Interferometer field 
of view) and to small scale structures down to 0.1 mas. 
As the visibility function is the Fourier transform of the 
image, a straightforward determination would require a 
huge number of image pixels (ss 10 4 x 10 4 ). 

As one of the goals of the present model is rapidity, 
we avoid this determination. Our image I(x, y) contains 
the central structure with about 10 3 x 10 3 pixels. All flux 
falling outside this image is integrated as a single value 
-Fout- Then, the non- normalised visibility V is given by 



V(u,v) 



1(0,0) +F out 
I(u/X,v/X) 



if u = v = 
otherwise 



(41) 



Fx = F x +F e x 



Ft 



(40a) 



where I is the Fourier transform of I. This equation states 
that the F out contribution cancels as soon as the baseline 
is non-zero, which means that the outer parts of the disc 
are fully resolved for any arbitrarily small baseline and 
would produce a visibility of if they were alone. Two ap- 
proximations are made: (i) the outer parts of the disc do 
not contain small-scale structures (ii) for too small base- 
lines Eq. H41fl is not valid. Our model does include neither 
irregularities nor steep variations in the outer parts of the 
disc so that approximation (i) is valid. The smallest base- 
line for which Eq. (|41|l is valid is determined by the field 
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Table 1. Parameters used for the fiducial model 



Parameters 


fiducial model 


r, (Re) 


2.0 


rmin (Re) 


6.0 


T-max (AU) 


100 


M* (M ) 


1.0 


M (Me/yr) 


l.OxlO" 8 


T* (K) 


4000 


in 





A v (mag) 


0.0 


,4 


0.2 


a 


lxlO" 2 



of view of the image 1 and we choose this field in order to 
produce correct results for the smallest baselines used in 
long-base optical interferometry (ss 10 m). 

4. Model analysis 

In this section study the influence of disc flaring hypothe- 
ses, compare our model with other authors' models, anal- 
yse the influence of different heating processes, show the 
influence of viscosity prescription and examine the influ- 
ence of some disc parameters on both the structure and the 
observables. For the sake of comparison, authors compute 
a fiducial model, typical of a T Tauri star; its parameters 
are displayed in Tabled 

From now on, we use model 2 in our study, unless 
specified otherwise (sect. |4~T1 and appendices). 

4.1. Importance of disc flaring hypotheses 

Figure |3| compares the structure of the fiducial T Tauri 
disc in three cases: the disc thickness is proportional to 
the scale height (model 1), the disc thickness is given by 
the reprocessing terms only (model 2), and the disc is flat 
(model 3). 

The major difference between models 1 and 2 is the 
presence of self-shadowing in the first one: at a radius of 
a few AUs, the surface temperature of model 1 drops be- 
cause the disc is not directly illuminated by the star (see 
middle panel of Fig.|3J). The cause is a decrease of the rel- 
ative scale height h/r (see right panel on the same figure), 
hence a similar behaviour of the disc surface because of 
the proportionality TL oc h. A contrario, model 2 assumes 
a thickness increasing faster than r (TC oc r 1+7 ), so the 
surface is always illuminated. 

Model 3 is a flat disc, therefore it catches much less 
light from the central star: the outer parts are much cooler 
(Figure |21 2], left panel). If the surface has the same tem- 
perature (cf. T c , middle panel), it produces a much smaller 
flux because the outer layer happens to be extremely thin. 



Such a disc has a lesser scale height because it is cooler 
(cf. h/r, right panel). 

As we shall see in Sect. 14.31 numerical simulations tend 
to pr ove that there is no sel f-shadowing effect in T Tauri 
discs. iD'Alessio et al.l l)l999j) establishes that even with a 
scale height flaring inwards, there is still enough material 
at large z to catch stellar light, so that all parts of the 
disc are illuminated. Moreover, we noticed that model 2, 
with the ad hoc flaring index, compares much better with 
other simulations than model 1. So, we shall continue our 
study with model 2. 

4.2. Viscosity prescription 

We compare the standard a viscosity prescription with 
the prescription d erived from laboratory experiments 



tne _jj prescription derived trom iabc 
bv lHure et ail (|2f)0llh is defined as 

vi = 0n K r 2 



(42) 



1 The visibility is the discretised Fourier transform of the 
image and the smallest non-zero baselines obtained are given 
by B = A/"!?, where •& is the field of the image. As a consequence 
Eq. 111! is only valid for baselines > B. 



and is a slowly varying factor ranging from 10~ 6 to 1CP 2 . 
We keep it constant over the disc. 

We chose = 10 -5 so that the viscosity v\ has the 
same order of magnitude as with the a-model in the cen- 
tral parts of the disc. This prescription gives much higher 
column density and disc mass in the outer parts of the 
disc, as displayed on Fig. Note that the impact of vis- 
cosity on the observables is small, since both viscous ef- 
fective temperature and reprocessing flux do not directly 
depend on it. In the mid-IR, the SED is only slightly af- 
fected, via a subtle modification of the disc thickness Tt, 
with an impact of the reprocessing temperature. Larger 
wavelength observations (> 100/zm) probe the outer parts 
of the disc, and should become more sensitive to the influ- 
ence of the viscosity law, because the flux emerging from 
these optically thin parts is proportional to kST. 

4.3. Comparison with other models 

We computed the two-layer structure of the fiducial model 
(Table and compare it model with others listed in 
Table [3 together with their main characteristics. FigureEl 
displays some physical conditions describing the radial 
structure of the fiducial disc, as forecast by these models. 
As a general result, despite of the approximation made, 
our model is in good agreement w ith previou s ones. 

Some discrepanc i es wit h the iBell et aP (|l997l^ and 
IChiang fc Goldreichl 1119971) models come from different 
heating hypotheses. IChiang Rr. Goldreichl (|l997h obtain 
lower values for temperatures in the inner parts of the 
disc, and higher mass columns since Sj/ is constant (see 
Fig. 15a. I & |5g.[ > . The reason is that they do not include 
viscous heating, predominant in the first AUs of T Tauri 
discs (see Fig. for r < 1 AU). In t he outer parts o f 
the disc, dominated by reprocessing, the lBell et al 1 lll997ft 
model differs significantly because it does not take this 
process into account. 

Our model and that of lD'Alessio et al"l 1 19991) present 
very close predictions in terms of central temperature, sur- 
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_ model 1 
. mode! 2 
model 3 




. model 1 
. model 2 
model 3 



a. 7] 



b. T c 



c. h/r 



Fig. 3. Influence of disc flaring hypotheses on the structure. Model 1: thickness proportional to the scale height, 
model 2: self-similar thickness, model 3: flat disc. Left: mid-plane temperature, middle: surface temperature, right: 
relative scale height. 




1.0 

r(u.a.) 




E 

g 10' 




a viscosity 

10-16 L p viscosity 



10~ 5 10 4 
Mm) 



a. Ti b. M d c. SED 

Fig. 4. Influen ce of the viscosity p rescription. Solid lines: a-prescription bv lShakura fc Sunvaevl l)l973[k dashed lines: (5- 
prescription bv lHure et alJ {2001). Right panel: mid-plane temperature vs. radius; middle panel: mid-plane temperature 
vs. radius; right panel: SED. 



face temperature or scale height. iHure k. Gallic |2001) 
already noticed that the flat vertically average disc model 
is a good approximate of an active disc; we have now 
demonstrated that a two-layer model is a faithful descrip- 
tion of discs both illuminated by the central star and 
heated by viscosity. However there are still some discrep- 
ancies: 



— We find important discrepancies in the region r < 
0.1 AU, because our opacities at intermediate tem- 
peratures (T > 1800 K) are different from those of 
iD'Alessio et all l|l999j) . 

These discrepancies do not affect much the results con- 
cerning SED and long-baseline visibility model fitting. 



— Our prediction of the disc thickness Ti is an under- 
estimation in the range 0.1-1 AU, so is the effective 
temperature of reprocessing T r . 

— Despite of a s maller effective temperature than in 
ID'Alessio et al. (1999). our SED is a slight overesti- 
mation in the range 20-100 /im. There is a likely rea- 
son: the two-layer structure tends to over-estimate the 
temperature of emission by the outer layer of the disc, 
hence the far-infrared excess. 



4.4. Location of heating processes 

Figure displays the radial profile of contributions of each 
heating process to the temperature of a typical T Tauri 
disc. The mid-plane temperature of the fiducial model is 
dominated by reprocessing for r > 2 AU and by viscous 
dissipation for r < 2 AU. However, in terms of effective 
temperature, the disc is dominated by reprocessing at any 
radial location, with a marginal contribution of viscosity 
for r < 0.2 AU. The outer layer is completely dominated 
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Table 2. Disc model characteristics compared in sec- 
tion 14.31 We mention whether stellar light reprocessing, 
viscous heating and convective transport are taken into 
account, whether the vertical transfer is solved or verti- 
cally averaged, and the viscosity prescription used. 



a b c d e 
reprocessing yes no yes yes no 

viscosity yes yes yes no yes 

convection no no no yes no 

vertical transfer no ID ID no ID 

viscosity prescription a, (3 a,/3 a — a 
a: present model 
b: paper II 

c: D'Alessio et al. f 19991 

H: Chiang Ooldreich (W97) 

c: Bell ct al. (1997) 



by stellar light reprocessing at any distance from the star. 
One important conclusion to draw is that the notion of 
a disc dominated by an heating process depends on what 
quantity we are interested in. 

4.5. Thermal flux vs. scattered light 

Simple disc models often ignore scattering when they fo- 
cus on the inner parts, but at large radius scattered light 
is dominant. A reason for this is that the thermal flux 
decreases exponentially with r, since it is proportional to 
B\(T(r)), while reprocessed flux is roughly in proportion- 
nal to r~ 2 — more exactly cx r~ 2 (f>. Figure[7]represents the 
contribution of thermal and scattered light to the emer- 
gent flux from a typical T Tauri disc, in different spectral 
bands. 

Let us first notice that the flux emerging from a stan- 
dard "one-layer" disc is dominated by scattered light for 
r > 0.15 AU in J, r > 0.5 AU in K, and r > 20 AU in 
N (contribution of the inner layer, Fig. \Q. Ignoring scat- 
tering in AU-scale predictions with the standard model 
proves out to be a bad approach in the near-IR. 

In a two-layer disc, the outer layer is hotter and pro- 
duces a larger thermal flux (contribution of the outer 
layer, same figure). Scattering is therefore dominant only 
at large radii: r > 2AU in J, r > 8AU in K, and 
r > 500 AU in N, so that ignoring it as a first approach 
can be an interesting simplification, when dealing with the 
AU-scale observations provided by IR interferometry, pro- 
vided that the field of view of the interferometer is limited 
(see Sect. 14.6. 1|) . In other words, the thermal emission is 
more extended in a two-layer disc, which implies for in- 
stance lower optical visibilities. 

4.6. Influence of disc parameters 

Figure |S] displays the influence of four parameters on the 
disc structure and the observables — the SED and the nor- 
malised visibility. 



The visibility curves are built using a point-like source, 
the star, with a flux i* 1 ,, over-imposed on an extended 
source -Fdisc- At large baselines, the visibility amplitude 
reaches a plateau yielding information on the flux ratio 
\V\ 2 w 1/(1 + F*/F di8c ) 2 . The higher the plateau, the 
smaller the ratio Fai sc /F*. The baseline B at which the 
plateau is reached gives the disc extent d w X/B. The 
higher the baseline B, the smaller the disc. 

At first sight the SED and visibilities present general 
features: the presence of a silicate feature in emission, due 
to the super-heated outer layer, and a steep drop of visi- 
bility when the wavelengths approaches the N band, also 
linked to the silicates: at 10 /im the outer layer becomes 
suddenly optically thicker and contributes more to the to- 
tal flux; since it is hotter and therefore presents a larger 
extent than the inner layer, the disc appears much larger. 

4.6.1. Albedo 

When the albedo varies, the outer parts of the structure 
and the SED are altered, because reprocessing is lowered. 
The effect on the structure of the inner parts is not visible, 
because reprocessing is not important there. The visible 
SED is only slightly changed because stellar light domi- 
nates scattered light by a factor > 10. 

While one usually ignores scattering in simple disc 
models, its incidence on the visibility curve can be im- 
portant. The scattered flux represents about AH/r of the 
stellar flux for a disc with a relative outer geometrical 
thickness H/r, i.e. w 5 to 10% in most discs. At near- 
infrared wavelengths and for moderate accretion rates, di- 
rect stellar flux and disc reprocessed flux are comparable, 
so scattered light cannot be ignored. Since reprocessing 
occurs at a large scale (like scattering, see Fig. , an ex- 
tended structure is resolved at very small baselines, hence 
the visibility drop near the origin. This effect might not be 
seen if the interferometer field of view is not large enough, 
hence a possible discrepancy up to 0.1 in normalised visi- 
bility. 

An important scattering also alters the variation of the 
visibility with wavelength by a differential effect: scatter- 
ing is much more important at short wavelengths, where 
it induces a larger drop in visibility. 

4.6.2. Inner disc radius 

As the inner truncation of the disc increases, the near- 
infrared part of the SED is depleted, because this part of 
the spectrum is mostly produced by the hot inner parts 
of the disc. 

The very long baseline (B > 400 m) visibility am- 
plitudes are higher, because the flux of the disc is lower 
compared to the flux of the star. The visibility plateau 
is reached for smaller baselines (ss 50-100 m) because 
the disc has a larger mean square dimension. For small 
baselines, the visibility is lower because the resolved outer 
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parts of the disc have a greater contribution compared to 
its still unresolved parts. 

The inner truncation also affects the variation of visi- 
bility with wavelength, because of a resolution effect. At 
small wavelengths, the resolution is high and visibilities 
follow the very long baseline trend and the visibility is 
higher for a large inner radius. At large wavelengths, the 
resolution is low and the visibility is lower for a large inner 
truncation. 

4.6.3. Viscosity 

As a increases, the viscosity becomes more efficient. The 
amount of material needed to produce a given accretion 
rate decreases. Therefore the optical thickness of the disc 
becomes lower and the mid-plane temperature follows the 
same trend. The discrepancy in the IR SED and visibil- 
ity amplitudes are small, because the effective tempera- 
ture of viscosity does not directly depend on a. However, 
when a increases, £ decreases, so that the thickness TL 
also decreases and less stellar light is caught by the disc. 
Therefore, the mid IR SED is depleted and the visibility is 
slightly higher. The SED is much more affected in the sub- 
millimetre, because it probes the optically thin outer parts 
of the disc, so that the flux is proportional to £ oc a -1 . 

4.6.4. Accretion rate 

As the accretion rate increases, the mid-plane tempera- 
ture is higher in the regions dominated by viscous heat- 
ing (r < 1-10 AU), but remains unchanged in the outer 
regions dominated by reprocessing. The near- and mid- 
infrared parts of the SED drastically change with accre- 
tion rate, since the viscous effective temperature is pro- 
portional to the accretion rate. At small accretion rates, 
the temperature inversion at the surface produces a sil- 
icate feature in emission at 10 fim; when M increases, 
the temperature inversion disappears, because the repro- 
cessing becomes secondary, and the feature disappears. 
The far-infrared SED does not change much for moderate 
accretion rates (10 _9 -10 -7 M Q /yr), because the regions 
that significantly contribute at large wavelength are dom- 
inated by stellar light reprocessing (almost) independent 
of the accretion. The slight decrease in the far infrared, 
when accretion increases, is due to the change of Tt with 
the amount of material. 

The visibility presents a lower plateau for larger accre- 
tion rates, because the flux from the disc becomes higher 
compared to the flux from the star. For large accretion 
rates (> 10~ 7 M Q /yr), the plateau is reached for smaller 
baselines because the mean angular size of the disc is 
larger (the hotter the disc, the larger the region of emission 
at a given wavelength). However at small baselines, discs 
with moderate accretion present a faster visibility drop 
because the resolved scattered flux is larger compared to 
the still unresolved inner parts. 



Table 3. Intcrfcromctric observations of FU Ori and 
T Tau stars 



star band instrument reference 

FU Ori H & K PTI/IOTA a 
T Tau N K PTI b 

SU Aur K PTI b 

a: Malbet et al. C1998) : Malbet fc Bcreer (2002) 
b: Akeson et al. f2000. 2002) 



Table 4. Parameters used for the models of our best fits. 
For T Tau North the SED and visibility fits are displayed. 





SU Aur 


FU Ori 


T Tau N 
(SED) 


T Tau N 
(V) 


r, (Rq) 


3.1 


4.0 a 


3.5 


3.8" 


r m in (Rq) 


8.2 


4.4 


3.6 


18.0 d 


r max (AU) 


200" 


150" 


80 a 


80" 


M* (M ) 


2.2 


1.0 


2.0 


2.0 d 


M (M /yr) 


2.0xl0 -7 


4.0 xlO -5 


1.4xl0" 7 c 


8.0xl0~ 7 


T, (K) 


5600 


5000 a 


4600 


4700 d 


in 


40 


40 







A v (mag) 


0.85 


1.4 


0.8 


0.8" 


A 


0.3 


0.2 


0.3 c 


0.2 d 


a 


2xl0- 46 


lxlO -1 b 


2xl0 -3 a 


lxl0" 2d 



a: parameter not constrained 

b: parameter loosely constrained 

c: (1 — A)M constrained 

d: parameters not separately constrained 



5. Comparison with observations 

We limit our study to TTS and FU Orionis stars for which 
both SEDs and optical visibilities are available. The sam- 
ple consists of three objects: FU Ori, SU Aur and T Tau 
North. Table 01 lists interferometric observations available 
for these stars. SEP measurements have been taken from 
iGezari et all l)l999|) . Table H lists the physical parameters 
of the best fits. 

Model- fitting of disc visib i lities h as already bee n car- 
ried out bv lMalbet fc Bergerl l|2002f) : lMalb"etl l|2002|k they 
consistently fitted the spectrum and K visibilities for 
FU Ori, with a self-similar flat disc model presenting an 
effective temperature oc r -0 ' 75 but were not able to repro- 
duce both H & K visibilities unless assuming a radial tem- 
perature law in r~ 0A to r~ 0,6 , that most physical accre- 
tion disc models cannot reproduce: in the case of FU Ori, 
dominated by viscous di ssipation, the expect ed temper- 
ature exponent is —0.75. lAkeson et alJ l|2002() separately 
fit SEDs and visibilities for T Tau and SU Aur but fail 
to find a set of parameters consistent both with interfero- 
metric data and SED. 

Figure displays observational data together with our 
best models fits. 

5.1. Fitting strategy 

We have identified various resulting characteristics of our 
modelled SEDs and visibilities that relate more specifically 



12 



Lachaume, Malbet & Monin: The vertical structure of T Tauri discs. Ill 



to physical parameters of the model. We describe hereafter 
how we use them to adjust our fits. 

— The level of the near-IR and mid-IR SED are very 
sensitive to the accretion rate (see Fig. l8d.|) . and more 
precisely to M^M (used to fit M*M). 

— The shape of near-IR and mid-IR SED is sensitive to 
the inner disc truncation (see Fig. 18b. \ (used to fit 

^"min ) • 

— The shape of the optical SED is sensitive to the stel- 
lar temperature, the stellar radius and the extinction 
(used to fit Ay, T* and r*). However, in massive discs, 
the shape of the optical SED is governed by the disc 
only and the values of r* and T* have little impact on 
the observables. 

— The shape of the far-IR SED is sensitive to the stellar 
temperature, the albedo and the parameter oj (used to 
fit uj and A). 

— The inclination of the disc implies a spread of visibility 
measurements at a given baseline modulus because the 
visibilities along major and minor axes then differ. We 
fit the inclination i using the spread in visibility. 

Visibilities have not been used to fit parameters di- 
rectly but were added as a constraint. The reason is 
that many parameters have an influence on visibilities, 
as shown in Fig. |H1 

We find values clo se to 0.2 for u> (0.1 in 
IChiang fc Ooldreichlll997|) . and albedo values close to 0.2. 
These parameters will no longer be free when we have 
included wavelength-dependent opacities. 



Moreover the low value for a baseline of 30 m a lso backs 
the second hypothesis. iMalbet fc Bergerl l|2002|) derived 
that the visibility measurements are compatible with a 
standard disc if a faint punctual source is located at 35 mas 
from the star. It might be interpreted as a hot spot in the 
disc. 

We are able to reproduce the SED until 20 /xm and 
the average of H and K measurements. The excess in 
the far IR is a well-know feature likely due to some phe- 
nomenon oco i ming_jn_the outer parts of the discs. For 
instance, lLodato fc BertirJ (|200lt) explain this excess with 
self-gravitation. 

Our model does not explain why the visibility value 
in the H-band at 100 m is much higher than in the K- 
band while the visibility curve in H is only slightly above 
the K curve. However our fit is consistent because of large 
uncertainties on visibility values. We indeed need more ac- 
curate measurements to disentangle betwe e n unc e rtainties 
and m odel inaccuracy. IMalbet fc Bergerl l|2002|) : IMalbet! 
(2002]) find that this higher value in H is consistent with 
an effective temperature distribution T v ~ r -°- 4 to -0.6 ^ 
that our model cannot reproduce. Figure l9b~1 displays the 
effective temperature of the disc: its law is T v ~ r -o.75^ m 
agreement with the active disc model. 

As expected, the structure of FU Ori (see mid-plane 
temperature curves Fig. 19b. I right panel) is dominated by 
viscous dissipation. Its effective temperature is also dom- 
inated by this process within a few tens of AUs from the 
star. 



5.2. SU Aur 

We were able to find a set of parameters to adjust simul- 
taneously the visibilities and the spectrum of this object 
(see Table ij. 

As expected with an accretion rate of 2 x 10 _7 M Q /yr, 
the structure of SU Aur is dominated by viscous heating 
within 10 AU from the star. (See the mid-plane tempera- 
ture curves on the right panel of Fig. 19a. Il On the contrary, 
the effective temperature of the disc is never dominated 
by viscous dissipation, as shown by effective temperature 
curves on the same figure. Therefore, the inner parts of 
the disc can be considered as an active disc in terms of 
structure but not in terms of observables. The reason for 
this contradictory behaviour is the following: T v « T r , but 
since T; ~ T v Ti + T r /2 with Tj 3> 1, we indeed have an in- 
ner temperature TJ w T v t u that is dominated by viscous 
heating. 

5.3. FU Ori 

The visibility measurements present a very large scatter. 
It may be produced by: (i) a large inclination of the disc 
(ii) a fully resolved structure. The first hypothesis can be 
discarded with scattered light images that do not show 
any jet, suggesting that the disc is seen almost pole on. 



5.4. T Tau North 

The available interferometric data for T Tau North in- 
clude the T Tau S outh component. We have used the 
lAkeson et all l|2002f) corrected data to adjust our model 
but we were unable to derive a unique set of physical pa- 
rameters to reproduce simultaneously the visibility and 
the SED of T Tau North. Therefore we present separate 
fits for the SED and the visibility, respectively displayed 
on Figs. El & IM] The best SED fit is obtained with a 
small inner disc truncation and a moderate accretion rate 
of 2.0 x 10~ 7 AfQ/yr. This fit overestimates the visibility 
by at least 2a. Visibility fits imply a large inner disc trun- 
cation and a larger accretion rate, inconsistent with the 
spectral data. 

We cannot exclude caveats in the visibility correction 
mentioned above. Hence we conclude that the physical 
parameters derived from the SED fit are the actual T Tau 
North's ones. Future VLTI or Keck Interferometer (KI) 
observations with a field of view of 50 mas should remove 
the T Tau N/S ambiguity. 

Like SU Aur, T Tau disc structure is dominated by 
viscosity within 10 AU from the star (see Fig. 19c. \ , but in 
terms of effective temperature both viscosity and repro- 
cessing must be taken into account. 
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6. Conclusion 

We have developed an analytical model for T Tauri accre- 
tion discs based on a two-layer approximation, and includ- 
ing the main heating processes: viscous dissipation, repro- 
cessing, and thcrmalization with the surrounding medium. 
The outer layer is directly heated by visible stellar light 
and the thermal flux from the inner layer, whereas the 
inner layer is heated by viscous dissipation and light re- 
processed by the outer layer. The strength of the model 
is an analytical prediction of the mid-plane temperature 
and, yet with less accuracy, of the flaring of the disc; it is 
a suitable tool in grasping the physical conditions taking 
place in these discs: it allows to predict easily the rela- 
tive importance of heating processes, the contribution of 
scattering, the influence of the viscosity prescription, etc. 

Despite of simplifications made in order to keep the 
model analytical, it compares well with other disc mod- 
els available in the literature. Its predictions in terms 
of mid-plane temperature, density scale height or disc 
mass are consistent w ith th ose of numerical models by 
iD'Alessio etall (|l998|h or bv lBell et all <|l997|) in the in- 
ner regions where this model is valid. 

For the first time, we are able to consistently fit the 
infrared spectra and the optical visibility of two young 
stars, SU Aur and FU Ori. A third object has been con- 
sidered (T Tau North), but we could not find a set of 
parameters that would fit simultaneously the SED and 
the visibility. This might result from the peculiar struc- 
ture of this object which is a triple system, and demon- 
strates that even a single interferometric measurement at 
one infrared wavelength can be a very strong constraint to 
interpret disc models. We therefore expect a breakthrough 
in disc physics understanding when new generation inter- 
ferometers, like the VLTI or the KI, are able to observe 
hundreds of young stellar objects. Providing an analytical 
description and a fast computation, while including essen- 
tial physical phenomena taking place in discs, our model 
will be a useful tool to interpret these forthcoming obser- 
vations. 

Concerning the influence of the viscosity on the model 
output, our results show that current instruments can- 
not significantly probe the influence of the viscosity law, 
because it has no direct impact on the flux emerging 
from optically thick regions. The dependency on viscos- 
ity observed in our model remains small and results from 
the variation of the disc geometrical thickness when the 
amount of material changes. It is only at larger wave- 
lengths, where the central regions are optically thin, that 
the mid-plane temperature can be probed, thus indirectly 
the viscosity The future Atacama Large Millimetre Array 
(ALMA) should be able to measure the column density 
and the mid-plane temperature of discs within 10 AUs 
from the star, and probe the influence of the viscosity. 

In forthcoming work, we will address the following 
points. 

— Including direct heating of the inner rim of the disc by 
the star will allow us to make predictions for more 



Table A.l. Five optical thickness regimes of the disc. 



regime 


TR 


TP 


T P 


Tp 


1 


> 1 


> 1 


» 1 


^> COS (f> 


2 


< 1 


> 1 


> 1 


<C COS (/> 


3 


< 1 


< 1 


> 1 


<C COS (/> 


4 


< 1 


< 1 


< 1 


<C cos cj> 


5 


< 1 


< 1 


< 1 


<§C cos 4> 



massive stars (Herbig Ae/Be systems) as shown by 
iDullemond et all ll200llh 

— The treatment of self-shadowing can also be improved: 
our model 1 assumes that the regions where the disc 
flares inwards are simply not illuminated at all, while 
vertical structure models tend to prove that there is 
still enough material at large heights to catch stel- 
lar radiation. Model 2 is neither satisfactory because 
the flaring index is taken ad hoc, though it gives quite 
faithful predictions. 

— A precise determination of the chemistry species in 
the outer parts of the disc implies to model the tem- 
perature profile in the outer layer ins t ead o f assum- 
ing an isothermal one. lAik awa et a"fl ll2002j) showed 
that the model by IChiang fc Goldreichl l)l997j) can- 
not reproduce observed element abundances in the 
ou ter part of the disc w hereas the vertical structure 
bv ID'Alessio et alJ l|l998|) can. 
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Appendix A: General n-disc solution 

In simplified disc models, dealing with optical thickness is 
always cumbersome, for five different mean optical thick- 
nesses of the inner layer have to be considered: Rosseland 
optical thickness tr, Planck thickness Tp, Planck-like 
thickness for the radiation of the outer layer Tp and the 
Planck- like thickness for the radiation of the star Tp. It 
leads to considering five different regimes summarised in 
Table IA.1I Most T Tauri discs present the regimes 1-4, 
while the fifth one is seldom encountered. For each of these 
regimes the mid-plane temperature follows different laws, 
expressed below: 



i^ + T^TR + i; 4 , 


(A.la) 


1 v a ' 


(A.lb) 


1 1 T 4 
i_L T 4 + iv_ T A 

2 2t p r 2t p a ' 




(A.lc) 


2t p r + 2t p + a ' 


(A.ld) 


1 Tp 2Tp ^,4 T v ^ 4 

2 Tp cos (f) 1 2t p 


(A.le) 



14 



Lachaume, Malbet & Monin: The vertical structure of T Tauri discs. Ill 

Derivating Eq. (|B.3|) in respect to r and introducing 

a k = (B.5a) 



where £ is a constant of the order of a few tenths related 
to the opacity local density exponent I (see Eq. I18fl . as 
defined in Eq. lf2"3|) . 

Some factors 2 in these expressions arise from our 
choice to deal with half-optical thickness (from mid-plane 
to outer layer). Instead of using four developments for the 
first four domains, we use an expression that follows the 
same asymptotic behaviours: 



T 



T t 4 F(2r|) , t4 



F(^)F(2r P ) 2F(2r P ) 
with 

F(x) = 1 — cxp(— x). 



(A.2) 



(A.3) 

In Eq. (|A~2)l . all F-terms can be interpreted using sim- 
ple radiative transfer arguments, except for 3 ^ that has 
been taken ad hoc. 

The different optical thicknesses of Eq. (|A.2() are de- 
termined from monochromatic opacities K\(p,T) by 



tr 



Skr(p, T, T c ) 
T,K P (p, T,T) 
Ek p (,o,T,T c ) 



(A.4) 
(A.5) 
(A.6) 



where the generalised Planck and Rosseland opacities of 
the material at (p, T) for black-body radiation at temper- 
ature T ra( j are given by: 



Kp 



Jo 



+°° dB x 
0T rac 



C^rad) dA 



lo 



+ 00 



dB x 



K X (p,T) 9T r! 



(T rad ) dA 



f °° K\{p, T)B\(T ia<i ) dA 



; +o ° 5 A (T rad )dA 



(A.7) 
(A.8) 



The opacity ratio in the outer layer (thermal radia- 
tion/stellar radiation) is then described as 



K P (p,T c ,T e ) 
K P (p,T e ,T*)' 



(A.9) 



Appendix B: Flaring index 

The mid-plane temperature is given by 
71 4 = 5>(r,T 5 ) 4 . 

k 

Using Eqs. l(TUll25)l one derives 
Ti = Ce(r) 2 r, 

where £ is a constant and e(r) = 
equation in to Eq. IjB.lfl we have 



(B.l) 



E (r) 8 = (Cr)- 4 £>(r,Cre(r) 3 ). 



(B.2) 

TL/r. Substituting this 
(B.3) 



The flaring is given by the variation of the relative 
geometrical thickness e(r) fEa. 126( 1: 



7 



r de 
e dr 



(B.4) 



dlogr ' 

4 
k 



dlogi 4 



dlogT; 
we derive 

y^7fct fc (r,Tj) 4 

k 

7 



5>(r,To 



with 



Ik 



4 + g k -b k 
8-2b k ' 



(B.5b) 



(B.6) 



(B.7) 



The flaring index of an optically thick a-disc is given 
in the main matter, by substituting the t^'s in Eq. IB.5| 
In the following we give the formulae in the case of an 
arbitrary optical thickness and of a /3-disc. 

B.l. a-disc of arbitrary optical thickness 

The expression for the temperature contributions t k are in 
the present case quite tangled (see Appendix lA"|l . though 
straightforward to derive, so we only give the results for 
the flaring index. We note 



xexp(— x) 
1 — exp(— x) 



(B.6 



the logarithmic derivative of F. The flaring associated to 
each heating process are: 



71 = 

72 = 

73 = 

74 = 

with 



V(l,m) (rK2r P ) - v (a£r) 


)+(2 + g(r)) 


Q(l,m) | 


[v(2tp) - i) 1 




) + 16 



V(l,m) (r ? (2rp)-7 ? (2rg))+4 
Q{l,m) ( ?7 (2Tp)-r ? (2T^)) + 14' 
V(l,m) (7 ? (2rp)-r / (2r|))+2 



Q(l,m) ( v (2tp)-ti(2t$)) 

1 

2 



16' 



V(l,m) = 2m + 31 + 1 
Q{l,m) = 4m - 61 - 4. 



(B.9a) 

(B.9b) 
(B.9c) 
(B.9d) 



(B.lOa) 
(B.lOb) 



B.2. Optically thick /3-d isc 

The switch to /3-prescription changes terms that involve 
column density: 



f(r)M 
3tt/3S1 



r- 



So ( 
r \ ir 



(B.H) 
(B.12) 
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Then, only the temperature contribution and flaring 
index related to viscous dissipation are changed: 



tiir) =^k t i / T g > p EoT v Ti 

t\\r) cx j(r) r ' T, , 
1-31 -2m 2 + 1 

7i 



16 + 21 -Am 16 + 2/- 4m' 
8.3. j3-disc of arbitrary optical thickness 



(B.13) 
(B.14) 
(B.15) 



The flaring indices are given by the same set of equations 
as in Appendix IB . II except for a change of Q(l,m): 

(B.16a) 
(B.16b) 



V(l, m) = 2m + 31 + 1 - (1 + l)g(r), 
Q(l,m) = 4m -21. 



Appendix C: Disc thickness on the outer edge 

In the main matter, we used an approximate g z cx z in 
the determination of the structure of the inner layer, be- 
cause most of the matter is concentrated within one or two 
scale heights, with hi <C r. In determining the vertical lo- 
cation of the disc surface with the same approximation, 
we noticed a large underestimation of H. The reason is 
that TL m 4-Qhi is not so small compared to r, so that the 
departure of g z from the linear law cannot be ignored. So, 
we are bound to considering the exact expression for the 
gravity field. 

In the outer layer, the Planck optical thickness r e , the 
surface density S e , and the opacity linked by 

T c = K C S . (C.l) 

If we use an isothermal layer, E can be substituted by 
p e h e where h e is the scale height of the outer layer and 
p c the density at its base. r c can be expressed as fcw, 
because we expect the incidence 4> of stellar radiation to 
be dominated by the flaring term. Therefore 



(C.2) 



To determine p c , we solve the vertical hydrostatic equilib- 
rium of the inner layer in the gravity field of the central 
star 



( \ ZC i 

9z{z) = -2hl 



which leads to 
p[z) = pi exp ^ 

with a scale height 
h(z) cx \ T{z) 1 







1 + 









-3/2 



1 - 



-1/2* 



(C.3) 



, (C.4) 



\ 



r 



-3/2 



(C.5) 



This allows us to write 

91 

2 



Pc = Pi exp ( — — S p . 



h r . 



(C.6) 



(C.7) 



S p and Sh being corrective terms due to the non-linearity 
of g(z). They are given in Sect. [TO (Eqs. fc |37|). 
Therefore 



Q 



2 . / KePihi T c 



(C.8) 



Equation (|35|) in Sect. 12.5.31 is then obtained after the 
substitution cj> = jQhi/r. 
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g. M d h. Q i. SED 

Fig. 5 . Co mparison between disc mod els: present work with a and (3 prescription. iD'Alessio et al.l ljl999t l. iBell et all 
l)l997|) and lChiang fc Goldreichl l|l997f) . The panels are: (a) mid-plane temperature (b) effective temperature associated 
to stellar irradiation (c) surface temperature (d) scale height (e) mean height of reprocessing (f) column density (g) 
cumulative disc mass (h) relative geometrical thickness of the disc (k) spectral energy distribution of the pole-on disc. 
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Fig. 6. Contribution of heating processes: radial profile of these contributions to the temperature of the disc (model 2). 
Left: mid-plane temperature, middle: surface temperature, right: surface temperature. 
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Fig. 7. Scattered light and thermal flux from a typical T Tauri disc (model 2) vs. the distance to the star. Left: J-band; 
middle: K-band; right: N-band. Solid lines: scattered light; dashed lines: thermal emission of the outer layer; dotted 
lines: thermal emission of the inner layer. 
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Fig. 8. Variation of the structure of a disc (first column), of its SED (second column), of its visibility curve (third 
column) , and of its visibility as a function of the wavelength (fourth column) when one of its parameters varies. 
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Fig. 9. Data and fits for YSOs (a: SU Aur, b: FU Ori, c,d: T Tau North's separate fits for SED and visibility). 
Left column: SED with contributions of the inner and outer layers and of the star; middle-left column: visibility 
amplitude vs. baseline, for the major and minor axes of the disc image; middle-right column: mid-plane temperature 
and contributions of reprocessing and viscous heating vs. radius; right column: effective temperature and contributions 
of reprocessing and viscous heating vs. radius. 



